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Abstract 

We report the observation of CP violation in the meson system. Using a novel technique 
for time-dependent measurements, we measure a non-zero value for the CP-violating ampli¬ 
tude sin2b at the 4.1 a level. We also report on precision measurements of the and 
lifetimes and the B^B^ mixing frequency Arrid obtained with the same technique, and on a 
first measurement of the time-dependent CP-violating amplitude in P° —> 7r’''7r“ decays. 
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1 Introduction 


CP violation has been a central concern of particle physics since its discovery in 1964 in the 
decays of mesons |]l|. An elegant explanation of the origin of CP violation was proposed 
by Kobayashi and Maskawa, as a complex phase in the three-generation CKM quark-mixing 
matrix P]. In this picture, measurements of CP-violating asymmetries in the time distributions 
of P° decays to charmonium hnal states are expected to be large and provide a direct test of 
the Standard Model of electroweak interactions [Q. 

We present measurements of time-dependent CP-asymmetries in samples of fully recon¬ 
structed B decays to charmonium-containing CP eigenstates {b —>■ ccs) and to the 7r’''7r“ final 
state. The data for these studies were recorded at the T(4S') resonance by the BABAR detector 
at the PEP-II asymmetric-energy collider at the Stanford Linear Accelerator Center. 

When the T(4S') decays, the P-wave B^B^ state evolves coherently until one of the mesons 
decays. In one of four time-order and flavor conhgurations, if the tagging meson Ptag decays 
first, and as a P°, the other meson must be a P° at that same time Aag- It then evolves 
independently, and can decay into a CP eigenstate Bcp at a later time tcp- The time between 
the two decays At = tcp — ttag is a signed quantity made measurable by producing the T(4S') 
with a boost /dy = 0.56 along the collision (z) axis, with nominal energies of 9.0 and 3.1 GeV 
for the electron and positron beams. The measured distance A 2 ; ftt'ycAt between the two 
decay vertices provides a good estimate of the corresponding time interval At; the average 
value of \ Az\ is /^ycrso 250/xm. 

We examine each of the events in the Bcp sample for evidence that the other neutral B 
meson decayed as a or a P° (flavor tag). The distribution f+(f-) of the decay rate when 
the tagging meson is a P°(P°) is given by 

p-AAP 

f±{At) = —--[1 ± Ssin{AmdAt) =f C cos(Amc/At)], (1) 

4r 


where At is the time between the two P decays, r is the B^ lifetime Am^ is the B^B^ 
mixing frequency 0 , and the lifetime difference between neutral P mass eigenstates is assumed 
to be negligible. The sine term in Eq. ^ is due to interference between direct decay and decay 
after mixing, and the cosine term is due to direct CP violation. The CP-violating parameters 
S and C are dehned in terms of a complex parameter A that depends on both P°P° mixing 
and on the amplitudes describing P° and P° decay to a common hnal state / 
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A difference between the P° and P° At distributions or a At asymmetry for either havor tag 
is evidence for CP violation. 

In the Standard Model A = for charmonium-containing b —> ccs decays, rjf is the 

CP eigenvalue of the state / and f3 = arg [ —KdKb/KdKb ] is an angle of the Unitarity Triangle 
of the three-generation CKM matrix |^. Thus, the time-dependent CP-violating asymmetry 

is 


Acp{At) = = -h/sm2^sin(AmB0 At). 


(3) 
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The analogous time-dependent CP-violating asymmetry in the decay —>■ 7T~^'n'~ arises 
from interference between mixing and decay amplitudes, and interference between the b —>■ 
uW~ (tree) and b ^ dg (penguin) decay amplitudes. If the decay proceeds purely through 
the tree process, the complex parameter A is directly related to CKM matrix elements and 
|A| = 1 and ImX = sin2a, where a = arg [-VtdV^l/VudV*,,]. 

However, recent theoretical estimates suggest that the contribution from the gluonic pen¬ 
guin amplitude can be signihcant 0 leading to |A| 7 ^ 1 and XmX = |A|sin2aeff, where tteg 
depends on the magnitudes and strong phases of the tree and penguin amplitudes. 

We also present precise measurements of the B^-B^ mixing frequency Am^ and the neutral 
and charged B lifetimes. These measurements use the same vertexing algorithm and At cal¬ 
culation as the measurements of the CP-asymmetries. In addition, for the Am^ measurement, 
the same flavor tagging algorithm as for the CP analyses is used. These measurements are 
amongst the most precise available and provide a good validation of the novel technique to 
study time-dependent B decays. 

In all analyses the values of the parameters under study (P lifetimes, Am^, sin2/3, and 
Ctttt ) were hidden to eliminate possible experimenter’s bias until event selection. At recon¬ 
struction method, and htting procedures were dualized and systematic errors were determined. 

2 The BaBAR Detector 

A detailed description of the BABAR detector can be found in Ref. [0. Charged particles are 
detected and their momenta measured by a combination of a silicon vertex tracker (SVT) 
consisting of hve double-sided layers and a central drift chamber (DCH), in a 1.5-T solenoidal 
held. The average vertex resolution in the z direction is 70 pm for a fully reconstructed 
B meson. We identify leptons and hadrons with measurements from all detector systems, 
including the energy loss (dP/dx) in the DCH and SVT. Electrons and photons are identihed 
by a Csl electromagnetic calorimeter (EMC). Muons are identihed in the instrumented hux 
return (IFR). A Cherenkov ring imaging detector (DIRC) covering the central region, together 
with the dP/da; information, provides K-tt separation of at least three standard deviations for 
B decay products with momentum greater than 250MeV/c in the laboratory. 


3 Measurement of B Lifetimes and Arud 

The measurements of the charged and neutral B lifetimes and the P°-P° mixing frequency 
Anid are based on a sample of approximately 23 million BB pairs. 

3.1 Exclusive B Reconstruction 

Samples of P° and B^ mesons Prec are reconstructed in the modes P° —> P+, 

P*^*^“a+, J/ipK*^ and P+ —» , 'i/;(2S')iC+. Charged and neutral D* candidates 

are formed by combining a P° with a n~ or 7 r°. P° candidates are reconstructed in the decay 
channels A+ 7 r“, A+tt^tt^, P+tt+tt^tt” and and D~ candidates in the decay channels 

and PsVr”. We reconstruct J/^lJ and in the decays to e+e“ and p+p“ and the 

'ip{2S) decay to J/'ip 7 r+ tt~. 
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Continuum e~^e~ —> qq background is suppressed by requirements on the normalized second 
Fox-Wolfram moment [§] for the event and on the angle between the thrust axes of i?rec and of 
the other B =-Bopp in the event. B candidates are identihed by the difference AE between the 
reconstructed B energy and the beam energy y/s/2 in the T(4S') frame, and the beam-energy 
substituted mass tjies calculated from a/s/2 and the reconstructed B momentum. We require 
rriEs > 5.2 GeV/c^ and | AE| < Sctae- The distributions of tues for selected B candidates in 30 
fb“^ is shown in Fig. |3.1| . 



Beam-Energy Substituted Mass (MeV/c 
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Figure 1: Beam energy substituted mass distribution for selected B^ (left) and (right) 
candidates. In 30 fb“^, we reconstruct 9400 B^ and 8500 B^ signal events. Average signal 
purities for uies > 5.27 GeV/c^, are 83 % and 85 % for B^ and respectively. 


3.2 At Measurement 

The decay time difference, At, between B decays is determined from the measured separation 
Az = Zrec — Zopp along the z axis between the i?rec and the i?opp vertices. This Az is converted 
into At using the T(4S') boost and correcting on an event-by-event basis for the direction of the 
B mesons. The resolution of the At measurement is dominated by the z resolution of the i?opp 
decay vertex. This vertex uses all tracks in the event except those incorporated in Brec- An 
additional constraint is provided by a calculated i^opp production point and three-momentum, 
determined from the three-momentum of the Brec candidate, its decay vertex, and the average 
position of the interaction point and the T(4S') boost. Reconstructed or A candidates are 
used as input to the fit in place of their daughters in order to reduce bias due to long-lived 
particles. Tracks with a large contribution to the are iteratively removed from the £t, until 
all remaining tracks have a reasonable fit probability. Gandidates with |Az| < 3.0 mm and 
(Taz < 400 /im are retained. For the measurement of the B lifetimes, we require that at least 
two tracks are included in the Ropp vertex . 

Two different parameterizations are used to model the decay-time difference resolution 
functions. In the measurements of Am^, sin2/3, and sin2Q;, the resolution function is approxi- 
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mated by a sum of three Gaussian distributions with different means 6k and widths cxfc, 
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where 6t is the difference between the measured and true At values. For the core and tail 
Gaussians, the widths ai ^2 = >S^i ,2 x (^At are the scaled event-by-event measurement error, 
a At, derived from the vertex hts. The third Gaussian, with a hxed width of = 8ps, 
accounts for less than 1% of outlier events with incorrectly reconstructed vertices. The three 
Gaussian resolution function is not suited for the measurement of the B lifetimes because 
strong correlations lead to increased statistical errors. Studies with Monte Garlo simulation 
and data show that the sum of a zero-mean Gaussian distribution and its convolution with an 
exponential provides a good trade-off between statistical and systematic uncertainties: 


'B{5t,aAt\a = {/i, s,k}) 





2s‘^a\J 
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The parameters a are the fraction h in the core Gaussian component, a scale factor s for the 
per-event errors a At, and the factor k in the effective time constant naAt of the exponential 
which accounts for charm decays. At outlier events are modeled the same way as in the 
three Gaussian resolution function. The resolution function parameters are assumed to be 
the same for all 5° and decay modes. This assumption is conhrmed by Monte Garlo 
simulation studies. The resolution functions differ only slightly between B^ and B^ mesons 
due to different mixtures oi D~ and mesons in the i?opp decays and we use a single set of 
resolution function parameters for both B^ and B^ in the lifetime hts. 


3.3 B Lifetime Results 

We extract the B^ and B^ lifetimes from an unbinned maximum likelihood ht to the At 
distributions of the selected B candidates. The probability for an event to be signal is esti¬ 
mated from uiEs fits (Fig. O) and the tues value of the i?rec candidate. In the likelihood, the 
probability density for the signal events is given by 

/ + 00 

g-Atl/r^(2r) TZ{At - At', aAt\a) d{At'), (6) 

■oo 

and the background At distribution for each B species is empirically modeled by the sum of a 
prompt component and a lifetime component convolved with the same resolution function, but 
with a separate set of parameters. The likelihood ht involves 17 free parameters in addition 
to the B^ and the B^ lifetimes: 12 to describe the background At distributions and 5 for 
the signal resolution function. The charged B lifetime tb+ is replaced with tb+ = r ■ r^o to 
estimate the statistical error on the ratio r = Tb+/tbo. 

We determine the B^ and B^ meson lifetimes and their ratio to be: 

tbo = 1.546 ± 0.032 (stat) ± 0.022 (syst) ps, 

tb+ = 1.673 ± 0.032 (stat) ± 0.023 (syst) ps, and 

tb+/tbo = 1.082 ± 0.026 (stat) ± 0.012 (syst). 
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These are the most precise measurements to date [Q and are consistent with the world 
averages The resolution function parameters are consistent with those found in a Monte 
Carlo simulation that includes detector alignment effects. Figure p.3| shows the results of the 
likelihood £t superimposed on the At distributions for and events. With the current 
data sample these measurements are still statistically limited. The dominant systematic errors 
arise from uncertainties in the description of the combinatorial background and of events with 
large At values, the use of a common time resolution function for B^ and B~^ and from limited 
Monte Carlo statistics. 



Ar(ps) At(ps) 


Figure 2: Left: At distribution for the B^ (top) and B~^ (bottom) events within 2 ct of the B 
mass with superimposed £t results. The single-hatched areas are the background contribu¬ 
tions and the cross-hatched areas represent the At outliers. The probability of obtaining a 
lower likelihood is 7.3%. Right: At distributions in data for the mixed and unmixed events 
(mEs(-Brec) > 5.27GeV/c^), with overlaid the projection of the likelihood £t (solid) and the 
background contributions (dashed). 


3.4 B Flavor Tagging 

The measurements of Am^, sin2/3 and sin2Q; require knowledge of the Ropp = -Stag flavor. We 
use the same tagging algorithm in the three analyses to determine the i?tag flavor from the 
charges of its decay products. 

The charge of energetic electrons and muons from semileptonic B decays, kaons, soft pions 
from D* decays, and high momentum charged particles is correlated with the flavor of the 
decaying b quark. Each event is assigned to one of four hierarchical, mutually exclusive tagging 
categories or has no flavor tag. A lepton tag requires an electron (muon) candidate with a 
center-of-mass momentum pcm > 1-0 (1.1) GeV/c. This efficiently selects primary leptons and 
reduces contamination due to oppositely-charged leptons from charm decays. Events meeting 
these criteria are assigned to the Lepton category unless the lepton charge and the net charge 
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of all kaon candidates indicate opposite tags. Events withont a lepton tag bnt with a non¬ 
zero net kaon charge are assigned to the Kaon category. All remaining events are passed to a 
nenral network algorithm whose main inpnts are the momentnm and charge of the track with 
the highest center-of-mass momentnm, and the ontpnts of secondary networks, trained with 
Monte Carlo samples to identify primary leptons, kaons, and soft pions. Based on the ontpnt 
of the nenral network algorithm, events are tagged as or and assigned to the NTl (more 
certain tags) or NT2 (less certain tags) category, or not tagged at all. The tagging power of 
the NTl and NT2 categories arises primarily from soft pions and from recovering nnidentified 
isolated primary electrons and mnons. The yields, efficiencies, pnrities and mistag rates w for 
each tagging category are listed in Table. |^. 

Table 1: Event yields for the different tagging categories obtained from fits to the mgs 
distribntions. The pnrity is qnoted for tties > 5.270 MeV/c^ and average mistag fractions Wi 
extracted for each tagging category i from the likelihood £t to the time distribution for the 
fully-reconstructed flavor eigenstate sample. 


Category 

Tagged 

Efficiency (%) 

Purity (%) 

w 

Lepton 

754 ± 28 

11.3 ±0.4 

97.1 ±0.6 

0.085 ±0.018 

Kaon 

2317 ±54 

34.8 ±0.6 

85.2 ±0.8 

0.167 ±0.014 

NTl 

556 ± 26 

8.3 ±0.3 

88.7 ± 1.5 

0.195 ±0.026 

NT2 

910 ±36 

13.7 ±0.4 

83.0 ± 1.3 

0.326 ±0.024 

Total 

4538 ± 75 

68.1 ±0.9 

86.7 ±0.5 



3.5 Mixing Result 


The value of Am^ is extracted from the tagged flavor-eigenstate B^ sample i?flav with a simul¬ 
taneous unbinned likelihood fit to the At distributions of both unmixed [B^B^) and mixed 
[B^B^ and B^B^) events. The PDFs for the unmixed (-I-) and mixed (—) signal events for 
the tagging category are given by 

7-f±(At| Arurf, Wj, hi) = —^-[1 ± (1 — 2wi) cos Am^^Af] 0 77(5t|aj), (7) 

Some resolution function parameters are allowed to differ for each tagging category to account 
for shifts due to inclusion of charm decay products in the tag vertex. The PDFs are extended 
to include background terms, different for each tagging category. The probability that a B^ 
candidate is a signal event is determined from a fit to the observed tues distribution for its 
tagging category. The At distributions of the combinatorial background are described with 
a zero lifetime component and a non-oscillatory component with non-zero lifetime. Separate 
resolution function parameters are used for signal and background to minimize correlations. 

The likelihood fit involves a total of 34 parameters, including Am^ (1), the mistag rates Wi 
and mistag differences AtCj = Wi{B^) —Wi{B^) (8), At resolution function parameters (9) and 
background parameters (16). We display the result of the likelihood fit by using the mixing 
asymmetry. 


*^mix {At) 
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unmixed 


A, 


unmixed 


{At) - NmixedjAt) 
{At) + Nmixed{At) 


( 8 ) 
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Figure 3: The asymmetry Amix{^t) between unmixed and mixed events as a function of |Af| 
overlaid with the result from the likelihood £t. 

If flavor tagging and At determination were perfect, the asymmetry as a function of At 
would be a cosine with unit amplitude. The amplitude is diluted by mistag probabilities and 
the At resolution. The At distributions of mixed and unmixed events, and their asymmetry, 
Amixi are shown in Fig. and along with projections of the £t result. The probability 
to obtain a smaller likelihood is 28 %. 

Systematic uncertainties in the Am^ measurement arise from various sources. The conver¬ 
sion of Az to At introduces an uncertainty (±0.007 ps“^) due to the limited knowledge of the 
PEP-II boost, the z length scale of BABAR and the i?rec momentum vector in the T(4S') frame. 
Systematic uncertainties related to the resolution function (±0.005 ps“^), are attributed to 
the choice of the parameterization, the description of outliers, and the capability of the resolu¬ 
tion model to deal with various plausible misalignment scenarios applied to the Monte Carlo 
simulation. The parameters of the background At distribution are left free in the likelihood 
£t, but systematic errors (±0.005 ps“^), are introduced by the uncertainty in signal proba¬ 
bilities, parameterization of the background At distributions and resolution function, and the 
small amount of correlated B^ background. Finally, statistical limitations of Monte Carlo 
validation tests (±0.004 ps“^), the full size of a (negative) correction obtained from Monte 
Carlo (±0.009 ps“^), and the variation of the B^ lifetime [Q (±0.006 ps“^) contribute. These 
contributions added in quadrature yield a total systematical error of 0.016 ps“^. 

In conclusion, the B^-B^ mixing frequency Airtd is determined with a new time-dependent 
technique from tagged fully-reconstructed B^ decays to be 

Arud = 0.519 ± 0.020(stat) ± 0.016(s|/st) ps“^. 

This preliminary result is one of the single most precise measurements available, and is con¬ 
sistent with the current world average and a recent BABAR measurement with a dilepton 


















sample |]OT. The error on Am^ is still dominated by statistics, leaving substantial room for 


further improvement as more data are accumulated. 


4 Measurement of sin2/? 

In 32 million BB, we extract sin2/? from a sample of fully reconstructed B^ decays (i?cp) to 
final states with rjf = —1 XciK^)i Vf = +1 ///(effective) = 

0.65 ± 0.07 (J/V'77*0 0. 

The sin2/9 measurement is made with a simultaneous unbinned likelihood £t to the At 
distributions of the tagged Bcp and Sgav samples. The At distribution of the former is 
given by Eq. |^, with |A| = 1. The i?flav sample evolves according to the PDF for B^ flavor 
oscillations as described in the previous section. The amplitudes for the Bcp CP-asymmetries 
and for the Pflav flavor oscillations are reduced by the same factor (1 — 2w) due to wrong tags. 
Both distributions are convolved with a common At resolution function and backgrounds are 
accounted for by adding terms to the likelihood, incorporated with different assumptions about 
their At evolution and convolved with a separate resolution function. Events are assigned 
signal and background probabilities based on the tues (all modes except J/'tpK^) or AE 
{J/'ipK^) distributions shown in Fig. |[ Separate At resolution functions parameters have 
been used for the data collected in 1999-2000 and 2001, due to the significant improvement in 
the SVT alignment. 

A total of 45 parameters are varied in the likelihood fit, including sin2/3 (1), mistag frac¬ 
tions w and differences Aw (8), parameters for the signal At resolution (16), and parameters 
for background time dependence (9), At resolution (3) and mistag fractions (8). The deter¬ 
mination of the mistag fractions and signal At resolution function is dominated by the large 
Pflav sample. The largest correlation between sin2/? and any linear combination of the other 
free parameters is only 0.13. We fix Tpo and Am^ [Q. 

Figure ^ shows the At distributions and Aqp as a function of At overlaid with the likelihood 
fit result for the rjf = —1 and rjf = +1 samples. The probability of obtaining a lower likelihood 
is 27%. The simultaneous fit to all CP decay modes and flavor decay modes yields 

sin2/3 = 0.59 ± 0.14 (stat) ± 0.05 (syst). 

The dominant sources of systematic error are the parameterization of the At resolution func¬ 
tion (0.03), due in part to residual uncertainties in SVT alignment, possible differences in 
the mistag fractions between the Bcp and Pfiav samples (0.03), and uncertainties in the level, 
composition, and CP asymmetry of the background in the selected CP events (0.02). The 
systematic errors from uncertainties in Am^o and r^o and from the parameterization of the 
background in the Pflav sample are small; an increase of 0.020 ps“^ in the value for Am^o 
decreases sin2/3 by 0.015. 

The large sample of reconstructed events allows a number of consistency checks, including 
separation of the data by decay mode, tagging category and Ptag flavor. The results of fits to 
some subsamples and to the samples of non-CP decay modes are shown in Table |^. For the 
latter, no statistically significant asymmetry is found. 

If |A| is allowed to float in the fit to the rjf = —1 sample, which has high purity and requires 
minimal assumptions on backgrounds, we obtain |A| = 0.93 ± 0.09 (stat) ± 0.03 (syst). The 
sources of systematic error are the same as in the sin2/9 analysis. 
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Figure 4: Left: a) distribution of tues for Bcp candidates having a in the hnal state; b) 
distribution of AE for J/xjjK^ candidates. Right: Number of rjf = —1 candidates {J/'ipKg, 
il){2S)Kl, and XciKf) in the signal region a) with a tag iV^o and b) with a tag 
Ngo, and c) the asymmetry [Npo — N^q)/{Nbo + N^o), as functions of At. The solid curves 
represent the result of the combined £t to all selected CP events; the shaded regions represent 
the background contributions. Figures d)-f) contain the corresponding information for the 
r]f = +1 mode 

Table 2: Number of tagged events, signal purity and observed CP asymmetries in the CP 
samples and control samples. Errors are statistical only. 


Sample 

Nta.g 

Purity (%) 

sm2p 

J/^Kl^{2S)KlxciK^s 

480 

96 

0.56±0.15 

J/^Kl 

273 

51 

0.70±0.34 


50 

74 

0.82±1.00 

Full CP sample 

803 

80 

0.59±0.14 

Rfiav non-CP sample 

7591 

86 

0.02±0.04 

Charged B non-CP sample 

6814 

86 

0.03±0.04 


The measurement of sin2/3 = 0.59 ±0.14 (stat) ±0.05 (syst) establishes CP violation in the 
B^ meson system at the 4.1cr level. This signihcance is computed from the sum in quadrature 
of the statistical and additive systematic errors. The probability of obtaining this value or 
higher in the absence of CP violation is less than 3 x 10“®. The corresponding probability for 
the rjf = —1 modes alone is 2 x 10“'^. 
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5 Measurement of sin2a 


We reconstruct neutral B mesons decaying to h^h'~, where h and h' refer to vr or in a 
sample of 33 million BB. The data set includes 30.4 fb“^ collected on the T(4S') resonance 
and 3.3 fb“^ collected below the BB threshold used for continuum background studies. 

20 
15 

^ 10 

i 5 
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I 80 
“60 
40 
20 


wigg (GeV/c^) 

Figure 5: Left: Distributions of tues and AE (unshaded histograms) for events enhanced in 
signal (a), (b) ttyt and (c), (d) iFvr decays. Solid curves represent projections of the maximum 
likelihood £t result, while dashed curves represent qq and ttyt Ktt cross-feed background. 
Shaded histograms show the subset of events that are tagged. Right: Distributions of At for 
events enhanced in signal tttt decays. Figures (a) and (b) show events (points with errors) 
with Rtag = B^ or B^. Solid curves represent projections of the likelihood £t, dashed curves 
represent the sum of qq and Kti background, and the shaded region represents the contribution 
from signal tttt. Figure (c) shows AnniAt) for data (points with errors), as well as £t projections 
for signal and background (solid curve), and signal only (dashed curve). 

We select B h'^h'~ candiates in the region 5.2 < mEs < 5.3 GeV/c^ and \ AE\ < 0.15 GeV 
and apply requirements on track multiplicity and event topology. The total number of events 
satisfying these criteria is 9741. This sample contains 97% background, mostly due to random 
combinations of tracks produced in > qq events. We extract signal and background 

yields for 7r’''7r“, K~^7r~, and K^K~ decays, and the amplitudes of the tttt sine (F^tt) and cosine 
(Gtttt) oscillation terms simultaneously from an unbinned likelihood £t. We parameterize the 
Ktt component in terms of the total yield and the CP-violating charge asymmetry Aktz = 
{.Nk--k+ — Nk+-k-){.Nk--k+ + Nk+-k-)- Background parameters are determined from tues and 
AE sideband regions. 

Discrimination between signal and background is based on tues, AP, and a Fisher dis¬ 
criminant T constructed from the scalar sum of the GM momenta of tracks and pho¬ 
tons (excluding tracks from the Prec candidate) and between pions and kaon tracks on the 
Gherenkov angle measurement from the DIRG. The inclusion of events with no flavor tag 
improves the signal yield estimates and provides a larger sample for determining background 
shape parameters in the likelihood fit. 
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There are 18 free parameters in the ht. In addition to the CP-violating parameters STr-n-, 
Ctttt, and Aktt (3), the ht determines signal and background yields (6), the background Kn 
charge asymmetry (1), and parameters describing the background shapes in m-Es, AP, and T 
(8). We hx r and Am^ |^. The At PDF for signal 7r’''7r“ decays is given by Eq. |^, modihed 
to include mistags and convolved with the signal resolution function. The At PDF for signal 
Pvt events takes into account P°-P° mixing and P° — K^K~ decays are parameterized as 
an exponential convolved with the resolution function. 

Figure shows distributions of mss and AP for events enhanced in signal decays based 
on likelihood ratios and the At distributions and CP asymmetry ^ 7 r 7 r(At) = (Ago (At) — 
A^o(At))/(Aso(At)-|-A^o(At)) for tagged events enhanced in signal tttt decays (approximately 
24 7171, 22 qq, and 5 Ptt events satisfy this selection). 

In conclusion, in a sample of 33 million BB, we hnd tttt, 217 ± 18 Kti, and d.dl®'! 
KK events. These yields are consistent with the branching fractions reported in Ref. We 
measure the following CP parameters: 

= O.OdlgJg (stat) ± 0.11 (syst), Aktt = —0.07 ± 0.08 (stat) ± 0.02 (syst), 

Ct,t, = —0.25^q; 47 (stat) ± 0.14 (syst). 

The correlation between and Ctttt is —21%, while Ak-k is uncorrelated with S'tttt and C-„t,. 
Systematic errors on Ct,.„, and Ak-k arise primarily from uncertainties in PDF shapes, 
tagging efficiencies and dilutions, r, and Am^. 

6 Summary and Outlook 

We have observed CP violation in the neutral B system at the 4.1 a level using a sample 
of fully reconstructed P° decays to CP eigenstates. With the same novel technique of time- 
dependent measurements, we have determined the B^ and P° lifetimes and the P°P° mixing 
frequency Am^ with high precision. In addition, we have presented the hrst measurement of 
the time-dependent asymmetry in P° —> 7r’''7r“ decays. All results are limited by the data 
sample size and we expect improved measurements from the rapidly growing BABAR data 
sample in the near future especially for the CP violating asymmetries. 
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